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Q_i' Abstract. We present a new method for solving the two-body Bethe-Salpeter equation in 

Minkowski space. It is based on the Nakanishi integral representation of the Bethe-Salpeter 
D I amplitude and on subsequent projection of the equation on the light-front plane. The method 

7^ , ■ is valid for any kernel given by the irreducible Feynman graphs and for systems of spinless 

particles or fermions. The Bethe-Salpeter amplitudes in Minkowski space are obtained. The 
^ , electromagnetic form factors are computed and compared to the Euclidean results. 

. 1 Introduction 

■ Bethe-Salpeter (BS) equation for a relativistic bound system was initially formulated in the 
! Minkowski space [1]. It determines the binding energy and the BS amplitude. However, in prac- 
tice, finding the solution in Minkowski space is made difficult due its singular behaviour. The 
singularities are integrable, but the standard approaches for solving integral equation fail. 
^ . To overcome this difficulty, Wick [2] formulated the BS equation in the Euclidean space, by 

^ [ rotating the relative energy in the complex plane fco — )• iA;4. This "Wick rotation" led to a well 
defined integral equation which can be solved by standard methods. Most of practical applica- 
tions of the BS equation have been achieved using this technique [Sj H] and recent developments 
make its solution a trivial numerical task [5] . Another method - the variational approach in the 
configuration Euclidean space - was recently developed in [6]. 

The binding energy provided by the Euclidean BS equation is the same than the Minkowski 
one. However, the Euclidean BS amplitude does not allow to calculate some observables, e.g. 
electromagnetic form factors. The integral providing the form factors contains singularities which 
are different from those appearing in the BS equation and whose positions depend on the mo- 
mentum transfer. Their existence invalidates the Wick rotation in the form factor integral. In 
terms of the Euclidean amplitude, the form factor can be obtained only approximately, in the so 
called static approximation. The error rapidly increases with the momentum transfer. To avoid 
this problem, the knowledge of the BS amplitude in Minkowski space is mandatory. 

Thus, fifty years after its formulation, finding the BS solutions in the Minkowski space is 
still a field of active research. 
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Solutions of the Bethe-Salpeter equation in Minkowski space 



Some attempts have been recently made to obtain the Minkowski BS amplitudes. The ap- 
proach proposed in [Tj [8] is based on the Nakanishi integral representation [9l |10] of the am- 
plitudes and solutions have been found for the ladder scalar case [71 [11] as well as, under some 
simplifying ansatz, for the fermionic one [12]. Another approach [13] relies on a separable ap- 
proximation of the kernel which leads to analytic solutions. Recent applications to the np system 
can be found in |14j . 

In the paper [15] we have proposed a new method to find the BS amplitude in Minkowski 
space and applied it to the system of two spinless particles. In [16] this approach was general- 
ized to the two-fermion case. Like in the papers [71 111^112]. our approach is based on Nakanishi 
representation of the BS amplitude. The main difference between our approach and those fol- 
lowed in [3 dH [12] is the fact that, in addition to the Nakanishi representation, we use the 
light-front projection. This eliminates the singularities related to the BS Minkowski amplitudes. 
Our method is valid for any kernel given by the irreducible Feynman graphs. 

In this paper we give brief review of the approach [15] to find the BS amplitude in Minkowski 
space, of its results and applications. 

The plan of the paper is the following. In order to present the method more distinctly, 
we consider first the case of zero total angular momentum and spinless particles. In sect. [21 
corresponding equation for the Nakanishi weight function is given. In sect. [3l it is applied to 
the ladder kernel and in sect. [4|- to the cross ladder one. In sect. [5] the method is generalized 
to the two-fermion system. Sect. [6|is devoted to application to the electromagnetic form factor, 
where advantage of the Minkowski space solution manifests itself in full measure. Sect. [7] contains 
concluding remarks. 



2 Equation for the v^eight function 

For a bound state of total momentum p and in case of equal mass particles, the BS equation 
reads 

= [(i + + / |^^^^'' ^»^(^''^)' 

where is the BS amplitude, iK the interaction kernel, m the mass of the constituents and k 
their relative momentum. We will denote by M = the total mass of the bound state, and 
hy B = 2m — M its binding energy. 

Our approach consists of two steps. In the first one, the BS amplitude is expressed via the 
Nakanishi integral representation [9l [10] : 



dz J d-/'- 



[/c^ + p ■ k z' + \M'^ — m? — 7' + ie 



3 • 



Notice that in this representation, the dependence on the two scalar arguments fc^ and p-k of 
the BS amplitude is made explicit by the integrand denominator and that the weight Nakanishi 
function (7(7, z) is non-singular. By inserting the amplitude ([2]) into the BS equation one finds 
an integral equation, still singular, for g. 

In the second step, we apply to both sides of BS equation an integral transform - light-front 
projection [15j - which eliminates singularities of the BS amplitude. It consists in replacing 
A: — 7- + j3 where w is a light-cone four-vector = 0, and integrating over j3 in infinite 
limits. We obtain in this way, a non-singular equation for the non-singular g{^,z). After solving 
it and substituting the solution in eq. ([2]), the BS amplitude in Minkowski space can be easily 
computed. 
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This leads (see ref. [15] for the detail of calculations) to the following equation for the weight 
function 51(7, z): 

9^^'^'^'^^' - [' dz' V{j,z;^',z')g{j',z'), (3) 



with 



7' + 7 + z^m^ + (1 — 

This is just the eigenvalue equation of our method. It is equivalent to the initial BS equation 
([1]). The total mass M of the system appears on both sides of equation ^ and is contained in 
the parameter = m? — |M^. As calculations [71 [8] show, g['y,z) may be zero in an interval 

< 7 < 7o- The exact value where it differs from zero is determined by the equation ([3]) itself. 

The kernel V, appearing in the right-hand side of eq. ([3]), is related to the kernel iK from 
the BS equation by 

V( > ') = ^ r -il{k + (3uj,p)d(3 

^^'^'^'^^ TT [(2 + A: + /3^^)2 - m2 + ie] [(2 - A; - /3w)2 - m2 + ze] ' 

For simple kernels K{k,k',p) given by a Feynman graph the integral ([5]) over k' is calculated 
analytically. The integral over /3 is also calculated analytically and it is expressed via residues. 
The singularities in the BS equation are removed by the analytical integration over /?. Equation 
([3]) is valid for an arbitrary kernel iK, given by a Feynman graph. The particular cases of the 
ladder kernel and of the Wick-Cutkosky model [21 [17] are detailed in the next section and for the 
cross ladder kernel ~ in the paper [18]. Once (7(7, z) is known, the BS amplitude can be restored 
by eq. p. 

The variables (7,2;) are related to the standard light front (LF) variables as ^ = k'j_, z = 

1 — 2x. The LF wave function (see its definition in [19], for instance) can be easily obtained by 

'ip{k^,x) = ^= -2- (6) 

V47r Jo y + ^2 ^ ^2 _ ^.(i _ 3.)^j2 

Eq. ([3]) can be transformed, in principle, to the equation for the LF wave function ip{k±,x), 
though this requires inverting the kernel in the left-hand side of ([3]). The initial BS equation ([T|), 
projected on the LF plane, can be also approximately transformed to the LF equation: 

fc^ -I- 777,2 ^\ , , fn^ f , , , I , 9. (fk\dx' , , 

-^^-M2jV'(fc±,x) = -^y ^(fcl,xOl^LF(fci,x^fc^,x,M2)^-;^^^ (7) 

with the LF kernel Vlf given, for ladder exchange, in ref. |15j . 

It is worth noticing that the LF wave function ([6|) is different from the one obtained by 
solving the ladder LF equation (jT]), as it was done e.g. in ref. [20]. The physical reason lies in 
the fact that the iterations of the ladder BS kernel (Feynman graphs) and the ladder LF kernel 
(time-ordered graphs) generate different intermediate states. The LF kernel and its iterations 
contain in the intermediate state only one exchanged particle, whereas the iterations of the 
ladder Feynman kernel contain also, many-body states with increasing number of exchanged 
particles (stretched boxes). This leads to a difference in the binding energies, which is however 
small [20]. Formally, this difference arises because of the approximations which are made in 
deriving eq. ([7|) from However, for a kernel given by a finite set of irreducible graphs, both 
BS ([1|) and LF ([7|) equations are already approximate and it is not evident which of them is 
more "physical". The physically transparent interpretation of the LF wave function makes it 
often more attractive. 
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Solutions of the Bethe-Salpeter equation in Minkowski space 



3 Spinless particles. Ladder kernel 

As illustration, we give here the kernel V{'y,z;'y',z') of equation ([3]) for the ladder BS kernel, 
which reads: 

We substitute it in eq. ([5]), then substitute ([5]) in dH and calculate the integrals. The details of 
these calculations are given in ref. |15j and one obtains: 



V[J,Z,J,Z) i ^(^^_^.y^_^/)^ ii-l<Z<z'<l ^' 



Wij,z,j',z') = --'-— ^ / -,dv (10) 



where W has the form: 

with a = (7^/(167rm^) and 

Do = + m^z^ + z'^)k^ 
D = v{l-v){l- z' )-/ + {!- z)[{l-v)n^ + vj'] 

+ vm'^ [(1 -v)il- z')z^ + vz''^{l - z)] + VK^{1 -z){l- z) \\^z-v{z- z) ] 

In the particular case we are considering here, the angular integral over v can be performed 
analytically and W obtains the simple expression: 

W{-f,z;j,z) = — (11) 



1 



(^+-^-)(2M--b+-^-) , , . 
(i-M(i-M + 



6± = {hi ± ^bl-^boh^ , 6o = (1 - z)m', 

61 = 7 + 7'-(l-z)^^-7'z-7z' + (l-z') [2;^m2 + (l-z2)«;2j ^ 

62 = -^{\- z')-{z- z')[{l- z){l- z')k^ + {z + z' - zz')m^]. 

One can show |15j that in the case = 0, which constitutes the original Wick-Cutkosky 
model [21 [17], the solution has the form 3(7, ^) = ^(7) g{z). The 7-dependence disappears from 
the equation which, after that, exactly coincides with the Wick-Cutkosky equation [2| [3l [T7]. 

Equation ([3]) has been solved by using the method explained in Appendix A, i.e. by expanding 
the solution g on a bicubic cubic spline basis. By keeping e/j = 10~^ fixed and varying the 
grid parameters to ensure four digits accuracy, we obtain for ji = 0.15 and /x = 0.5 and unit 
constituent mass (m = 1) the values displayed in tabled) They correspond to 7max = 3, A''^ = 12, 
Nz = 10. With all shown digits, they are in full agreement with the results we have obtained, 
similarly to [2^, by using the Wick rotation and the method of [5]. Increasing e/j to 10~^ changes 
at most one unit in the last digit. This demonstrates the validity of our approach. 

The weight function g for a system with // = 0.5 and B = 1.0 is plotted in fig. [TJ It has been 
obtained with er = 10~^ and the same grid parameters than in table [TJ Its 7-dependence is not 
monotonous and has a nodal structure; the z-variation is also non trivial. We have remarked 
a strong dependence of g{^,z) relative to values of the er parameter smaller than ~ 10~^, in 
contrast to high stability of corresponding eigenvalues. However the corresponding BS amplitude 
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Table 1. Coupling constant values as a function of the binding energy for /i = 0.15 and fi — 0.5 obtained with 

7maa- ^3, N.y = 12, TV, = 10 and en. = 10"^ 



B 


a{n = 0.15) 


a(/i = 0.50) 


0.01 


0.5716 


1.440 


0.10 


1.437 


2.498 


0.20 


2.100 


3.251 


0.50 


3.611 


4.901 


1.00 


5.315 


6.712 




Figure 1. Nakanishi function g{'y,z) for fi = 0.5 and B = 1.0. On left - versus 7 for fixed values of z and on 
right - versus z for a few fixed values of 7. 



^ and LF wave function ip, obtained from g{'y, z) by the integrals ([2]) and ([6]), show the same 
strong stabihty as the eigenvalues. 

The BS amplitude in Minkowski space in the rest frame p = is shown in fig. [2j The k- 
dependence is rather smooth but the feo-dependence, due to poles of the propagators in ([T]), 
exhibits a singular behaviour at /cq = i (e/c it i.e. moving with k and M. 




Figure 2. BS amplitude <l>{ko,k) for fi = 0.5 and B = 1.0. On left versus k for fixed values of fco and on right 
versus fco for a few fixed values of fc. 

Note that our solution gives also the BS amplitude in Euclidean space, by substituting in ([2]) 
ko = iki. It is indistinguishable from the one obtained by a direct solution of the Wick-rotated 
BS equation [15]. 
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Solutions of the Bethe-Salpeter equation in Minkowski space 



4 Spinless particles. Cross ladder kernel 

Non-ladder effects, within the same model, using Feynman-Schwinger representation, were con- 
sidered in ref. [22]. In this work the full set of all irreducible cross-ladder graphs in a bound 
state calculation was included. Refs. [231 [24] contain results on the binding energy found by 
solving the BS equation for (L+CL) in Euclidean space. In the LFD framework, binding energy 
for (L+SB) kernel is calculated in [25], whereas the (L+CL+SB) contribution is incorporated 
in |24] . In [26] the effect of the cross-ladder graphs in the BS framework was estimated with 
the kernel represented through a dispersion relation. Non-ladder self-energy effects have been 
incorporated in [27] [28 ] fTT] . 

p" 

Pi \ ^ P> 

Pj - p" "■ ■■■ . Pi - P" 

p, — — p'2 

p;-p, + p" 

Figure 3. Feynman cross graph. 

As mentioned, the equation ^ is valid for any kernel. Derivation of V in this equation for the 
cross-ladder BS kernel, shown in fig. [S] is quite similar but more lengthy, since the kernel itself is 
more complicated. We have first to calculate to the kernel K^^^\k,k' ,p) corresponding to the 
diagram in fig. [S] substitute the result in ([5]), then in (jl]) and find in this way the cross-ladder 
contribution to the kernel 1^(7, 7', 2') in equation ([3]). The full kernel - including ladder and 
cross-ladder graphs - will be written in the form: 

1/(7, z; 7', z') = y(^) (7, z- 7', z') + (7, z; 7', z). 

The ladder kernel is given in refs. |15[ I19|. The cross-ladder contribution was calcu- 

lated in the paper [18] . 

We compare the results obtained in the BS approach with the equivalent ones found in Light- 
Front Dynamics (LFD). The latter incorporates the graphs with two bosons in fiight forming 
six cross boxes and two stretched boxes. They were also calculated in [18| . 

The binding energy B as a, function of the coupling constant a is shown in figures |4] for 
exchange masses /i = 0.15 and /U = 0.5 respectively and m = 1. 

We see that for the same kernel - ladder or (ladder + cross-ladder) - and exchange mass - 
= 0.15 or 11 = 0.5 ~ the binding energies obtained by BS and LFD approaches are very close 
to each other. The BS equation is slightly more attractive than LFD. At the same time, the 
results for ladder and (ladder -l-cross-ladder) kernels considerably differ from each other. The 
effect of the cross-ladder is strongly attractive. Though the stretched box graphs are included, 
its contribution to the binding energy is smaller than 2% and attractive as well. 

The zero binding limit of fig. |3] deserves some comments. It was found in [20] that for massive 
exchange, the relativistic (BS and LF) ladder results do not coincide with those provided by the 
Schrodinger equation and the corresponding non relativistic kernel (Yukawa potential) even at 
very small binding energies. Their differences increase with the exchanged mass ^ and do not 
vanish in the limit I? — t- 0. We have displayed in fig. [5] a zoom of fig. U] (right) for small values 
of B. The cross ladder and stretched box diagrams reduce the differences but are not enough to 
cancel it. 

We see that the cross-ladder contribution, relative to the ladder one, results in a strong 
attractive effect. The BS and LFD approaches give very close results for any kernel, with BS 
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Figure 4. On left: Binding energy B vs. coupling constant a for BS and LF equations with the ladder (L) kernels 
only and with the ladder +cross-ladder (L+CL) one for exchange mass /i — 0.15. On right: The same as on left 
but for exchange mass fx = 0.5 and, in addition, binding energy B for LF equation with the ladder +cross-ladder 
+stretched box (L+CL+SB) kernel. 



equation being always more attractive. These approaches differ from each other by the stretched- 
box diagrams with higher numbers of intermediate mesons. Our results indicate that the higher 
order stretched box contributions are small. This agrees with direct calculations in LFD of 
stretched box kernel with two-meson states [29] and with calculations of the higher Fock sector 
contributions [30] in the Wick-Cutkosky model. Calculation in LFD of binding energy with the 
stretched box contribution (L+CL+SB) and its comparison with (L+CL) also shows that the 
stretched box contribution is attractive but small. 



0.5 





LF Ladder 

LF Ladder+CL+SB 

Schrodtnger (Yukawa) 








(x=0.50 



0.01 0.02 0.03 

B 



Figure 5. Zoom of figure |4] (right) in the zero binding energy region. The ladder and (ladder +cross ladder 
+stretched box) results obtained with the LF equation are compared to the non relativistic ones (Schrodinger 
equation with Yukawa potential). 



The comparison of our results with those obtained in [22] , evaluating the binding energy Baii 
for the complete set of all irreducible diagrams, shows that the effect of the considered cross 
ladder graphs, though being very important, represent only a small part of the total correction. 
Thus for /i = 0.15 and a = 0.9 the corresponding binding energies obtained with BS equation 
are Bl ~ 0.035, ^l+cl ~ 0.06 and Baii ~ 0.225. 

Qualitatively, our large attractive effect of CL (for /i = 0.15 and /i = 0.5) is in agreement 
with |23j (where [i=V). For /u = 0.15, our numerical results for binding energy, both for the BS 
equation in Minkowski space and in the LFD approach, are close to those found in [21] (within 
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Solutions of the Bethe-Salpeter equation in Minkowski space 



accuracy of extracting the data from plots). Our results are smaller than the ones found in [26] 
by a factor 3. 

5 Two fermions 

We have considered the following fermion m) - meson fj) interaction Lagrangians: 

(i) Scalar coupling df^^ = g^cj)^ for which Fa = ig 

(ii) Pseudoscalar coupling c\^f = ig^^^cp^ for which Fa = —gjb- 

(ii) Massless vector exchange £^-^1 = ^jt^V^ ^ with Fa = igj'^ and TI^j/ = —ig^u/Q"^ as vector 
propagator. 

Each interaction vertex has been regularized with a vertex form factor F[k — k') by the 
replacement g — )■ gF{k — k') and we have chosen F in the form: 

F{q) = -J^^ -. (12) 

Let us first consider the case of the scalar coupling and the corresponding ladder kernel. The 
BS equation for the amplitude reads: 



<^{k,p) 



i{m + 2P + k) 



{i^p + kY — rn? + ie 



d'fc' ^ {-ig^)F\k-k') 



{2ttY ' {k-k>Y - ^j? +ie 



i{m — ^p + k) 



{^p — ky — m? + ie 

(13) 

In the case of J'" = 0"*" state, the BS amplitude has the following general form: 

${k,p) =Si(t>l+ S2(t>2 + S3<p3 + 54.^4 (14) 

where Si are independent spin structures (4x4 matrices) and (pi are scalar functions of /c^ and 
p-k. The choice of Si is to some extent arbitrary. To benefit from useful orthogonality properties 
we have taken 

Si =75, S2 = ^P15, ^3 = -^^75-^^75, Si = ^(J^^p^ky-i^, 

where cj^jy = §(7^7i/ — 71^7^). The antisymmetry of the amplitude (fHl) with respect to the permu- 
tation 102 implies for the scalar functions: 4'i,2,A{k,p) = (pi,2,4{—f^,p), 4''i{k,p) = —(j)3{—k,p). 
A decomposition similar to (I14p was used in [12] to solve the BS equation for a quark-antiquark 
system but the solution was approximated keeping only the first term Si(j)i. 

We substitute (fUj) in eq. (fT3l) . multiply it by Si and take traces. As we will see, the kernel 
in the resulted equation, in contrast to the spinless case, is still singular. These singularities are 
integrable numerically. They do not prevent from finding numerical solution, but they reduce 
its precision. This can be avoided by a proper regularization of equation, multiplying both sides 
of it by the factor 

'^^''^P' [(2 + ky - L2 + ie] [(I - fc)2 - L2 + ie] ^ ' 

This factor has the form of a product of two scalar propagators with mass L. It plays the 
role of form factor suppressing the high off-mass shell values of the constituent four-momenta 
ki2 = (2^^)^ ^^'^ tends to 1 when L — )■ 00. In this way, we get the following system of equations 
for the invariant functions (hi: 



ri{k,p) (l)i{k,p) 



ri{k,p) 



[(I + /c)2 - ie][(| - kf - + ie] 
f d'^k' ig^F\k-k') ^ 
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Since r]{k,p) / 0, the equation thus obtained is strictly equivalent to one where i]{k,p) is 
cancelled. We will see however that, due to the presence of the t] factor, the LF projection 
modifies the resulting kernels which become less singular functions. 

The coefficients Cij are determined by traces. We do not give here their explicit form, which 
can be found in [16| . 

Then we represent each of the BS components (j)i{k,p) by means of the Nakanishi integral ([2]) 
and, similarly to the scalar case, apply the light-front projection to the set of coupled equations 
for the corresponding weight functions gi{'y,z). As explained in sect. El this projection, which is 
an essential ingredient of our works [151 [T8| I31j . consists in replacing k ^ k + in eq. ()16p 
and integrating over /3 in all the real domain. 

The technical details of the light-front projection are similar to those given in ref. [15]. We 
obtain in this way a set of coupled two-dimensional integral equations: 

poo rl poo pi 

H j^^dz'V3{^,z-i,z')gH,z) = Y,j^ j dz'V,^^{^,z-i,z')g,{i,z') (17) 

The kernel and also for all types of couplings and states are given in jl6j . These kernels 
depend on the parameter L. Closer is L to m, smoother is the kernel and more stable are the 
numerical solutions. However the weight functions gi('j,z) as well as binding energies provided 
by (fTT]) are independent of L. To avoid spurious singularities in (flTl) due to the r] factor (fT5]l . 

must be larger than ^ m^, what is fulfilled for L > m. In practical calculations we have 

taken m = 1 and L = 1.1. 

The kernel Vg is finite and it vanishes for z = ±1. For a fixed values of j,z and 7', is a 
continuous function of z' with a discontinuous derivative at z' = z. 




^).5 0.6 0.7 0.8 0.9 1 -0.005^ ^ 

z' z' 



Figure 6. Left: kernel matrix elements Vi4(7, 2; 7', z') for z = 0.7,0.9,0.95 as a function of z' and fixed values 
of 7,7'. The discontinuity is finite at a fixed value of z but diverges when 2 — > 1. Right: regularized kernel V14, 
plotted for the same arguments and L — 1.1m. 



As already mentioned, without ?7-factor, most of the kernel matrix elements V^j are singular. 
Namely, they are discontinuous at z' = z. In some cases - like e.g. V14 displayed in fig. [6] (left) - 
the value of the discontinuity, although being finite at fixed value of z, diverges when z ±1. 
This creates some numerical difficulties when computing the solutions 51(7, z) in the vicinity of 
z = ±1. They can be in principle solved by properly taking into account the particular type 
of divergence. However the latter may depend on the particular matrix element, on the type of 
coupling, the quantum number of the state and other details of the calculation. 

For T] 7^ 1, i.e., for a finite value of L in (jlSp . the ^'-dependence of the regularized kernels 
is much more smooth and therefore better adapted for obtaining accurate numerical solutions. 
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Solutions of the Bethe-Salpeter equation in Minkowski space 



In fig. [6] (right) we plotted the regularized kernel Vf^ as a function of z' for the same arguments 
7, z, 7' and parameters than in fig. [6] (left), where it was calculated without the r]{k,p) factor. As 
one can see, the kernel is now a continuous function of z'. A discontinuity of the first derivative 
however remains at z' = z. 

We would like to emphasize again that despite the fact that the non-regularized and regu- 
larized kernels are very different from each other (compare e.g. figs. [6] left and right) and that 
the regularized ones strongly depends on the value of L, they provide - up to numerical errors 
— the same binding energies and weight functions 54(7, z). We construct in this way a family of 
equivalent kernels. 



Table 2. Left: Coupling constant as a function of binding energy B for the J — Q state with scalar (S), 
pseudoscalar (PS) and massless vector (positronium) exchange kernels. The vertex form factor is yl = 2 and 
the parameter of the 77 factor L = 1.1. Right: Coupling constant as a function of binding energy B for the 
positronium J = Q state in BS equation in the region of stability without vertex form factor [A — !> 00), i.e. g < n. 
They are compared to the non relativistic results. 





s 


PS 


positronium 




0.15 


0.50 


0.15 


0.50 


0.0 


B 


9' 


9' 


9' 


9' 


9' 


0.01 


7.813 


25.23 


224.8 


422.3 


3.265 


0.02 


10.05 


29.49 


232.9 


430.1 


4.910 


0.03 


11.95 


33.01 


238.5 


435.8 


6.263 


0.04 


13.69 


36.19 


243.1 


440.4 


7.457 


0.05 


15.35 


39.19 


247.0 


444.3 


8.548 


0.10 


23.12 


52.82 


262.1 


459.9 


13.15 


0.20 


38.32 


78.25 


282.9 


480.7 


20.43 


0.30 


54.20 


103.8 


298.6 


497.4 


26.50 


0.40 


71.07 


130.7 


311.8 


515.2 


31.84 


0.50 


86.95 


157.4 


323.1 


525.9 


36.62 



B 


9nr 


9ls 


0.01 


2.51 


3.18 


0.02 


3.55 


4.65 


0.03 


4.35 


5.75 


0.04 


5.03 


6.64 


0.05 


5.62 


7.38 


0.06 


7.95 


8.02 


0.07 


11.24 


8.57 


0.08 


13.77 


9.06 


0.09 


15.90 


9.49 



The solutions of eq. (fT7|) have been obtained using the techniques detailed in Appendix A. 
We have computed the binding energies and BS amplitudes, for the J = 0^ two fermion system 
interacting with massive scalar (S) and pseudoscalar (PS) exchange kernels and for the fermion- 
antifermion system interacting with massless vector exchange in Feynman gauge. In the limit of 
an infinite vertex form factor parameter yl — t- 00, the later case would correspond to positronium 
with an arbitrary value of the coupling constant. All the results presented in this section are 
given in the constituent mass units (m = 1) and with L = 1.1. 

The binding energies obtained with the form factor parameter yl = 2 are given in left table 
[21 For the scalar and pseudoscalar cases, we present the results for = 0.15 and fi = 0.50 boson 
masses. They have been compared to those obtained in a previous calculation in Euclidean 
space [32] using a slightly different form factor, which differs from our one by a factor. Once 
taken into account this correction, our scalar results are in full agreement (four digits) with 
[32] . The pseudoscalar ones show small discrepancies (ss 0.5%). We have also computed the 
binding energies by directly solving the fermion BS equation the Euclidean space using a method 
independent of the one used in |32j . Our Euclidean results are in full agreement with those given 
in the table [2j 

The B{g'^) dependence for the scalar and pseudoscalar couplings is plotted in figs. [71 Notice 
the different scales of both dependences. The pseudoscalar binding energies are fast increas- 




ing functions of and thus more sensitive to the accuracy of numerical methods. This sharp 
behaviour was also exhibit when solving the corresponding LF equation |33j . 

It is worth noticing that the stability properties of the BS J'^ = O"*" solutions for the scalar 
coupling are very similar to the LF ones. In the latter case, we have shown [341135] the existence of 
a critical coupling constant Qc below which the system is stable without vertex form factor while 
for g > gc, the system "collapses", i.e. the spectrum is unbounded from below. The numerical 
value was found to be ac = 3.72 \34:\ 135] . which corresponds to gc = 6.84. Performing the same 
analysis than in our previous work ~ eq. (71) from \3A\ [35] - we found that for BS equation the 
critical coupling constant is gc = 27r, in agreement with [32]. The 10% difference between the 
numerical values of gc is apparently due to the different contents of the intermediate states in 
the two approaches. The ladder BS equation incorporates effectively the so-called stretch-boxes 
diagrams which are not taken into account in the ladder LF results. 

The positronium case deserves some comments. First we would like to notice that in our 
formalism, the singularity of the Coulomb-like kernels in terms of the momentum transfer 
l/{k - k'f is absent. This is a combined consequence of the Nakanishi transform ([2]) - which 
allows to integrate over k' analytically in the right hand side of the BS equation ()16p - and of the 
consecutive LF projection integral. After this integration, the Coulomb singularity does not any- 
more manifest itself in the kernel. This can be explicitly seen in the kernel of the Wick-Cutkosky 
model obtained in eq. (22) of our previous work |15] . 

A second remark concerns the A dependence of the positronium results. Using the methods 
developed in |341 [35] we found that in the BS approach with ladder kernel there also exists a 
critical value of the coupling constant gc = tt. Note that, as in the scalar coupling, the very 
existence and the value of this critical coupling constant is independent on the constituent (m) 
and exchange /i masses but depends on the quantum number of the state. 

The ground state positronium binding energies without vertex form factor are given in the 
left table [2] for values of the coupling below gc, Nonrelativistic results gffj^ = Siry/B/m are 
included for comparison. One can see that the relativistic effects are repulsive. 

These results are displayed in fig. [8] (left), black solid line, and compared to the binding 
energies obtained with two values of the form factor parameter A = 2 (dashed) and yl = 5 (dot- 
dashed). The stability region is limited by a vertical dotted line at g = gc = tt. Beyond this value 
the binding energy without form factor becomes infinite and we have found B{g — >■ gc) ~ 0.10. 

The inclusion of the form form factor has a repulsive effect, i. e. for a fixed value of the coupling 
constant it provides a binding energy of the system which is smaller than in the yl — >• oo limit (no 
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Positronium J=0 




50 



Figure 8. Left: Binding energy for J=0 positronium state versus (black solid line) in the stability region 
g < Qc — n. Dashed and dotted-dashed curves correspond to the results for increasing values of the vertex form 
factor parameter A. They are compared to the non relativistic results (red solid line). Right: yl-dependence of g^ 
for for J — positronium state for fixed values of binding energies B = 0.1 and B — 0.5. 



cut-off). This is also illustrated in fig. [8] (right) where we have plotted the A dependence of for 
two different energies. One can see that the value of the coupling constant to produce a bound 
state is a decreasing function of A. The size of the effect depends strongly on the binding energy 
but for both energies the asymptotics is reached at A ^ 20. This behaviour is understandable 
in terms of regularizing the short range singularity of — interactions. 





a \ 2 3 i .0.5 0.5 1 

Y z 
Figure 9. Left: Nakanishi weight functions v.s. 7 for z — 0.6, for scalar exchange for yl = 2, L = 1.1, /i = 0.15 

and n = 0.5. Right: Nakanishi weight functions v.s. z for 7 — 0.54. 

We present in figs. [9] some examples of the Nakanishi weigh functions gi{'y,z). They corre- 
spond to a i? = 0.1 state with the scalar coupling and the same parameters A = 2, fi = 0.50 
than in table [2j In the left figure is shown the 7-dependence for a fixed value of z and in the 
right figure the z-dependence for a fixed 7. Notice the regular behaviour of these functions as 
well as their well defined parity with respect to z - 51,2,4 are even and gs is odd. As in the scalar 
case, the e/j-dependence of gi is more important than for the binding energy. 

Corresponding BS amplitudes (pi are displayed in figs.[10l The figure [10] (left) represents the 
ko dependence of (pi for a fixed value of | k |= 0.2. They exhibit a singular behaviour which 
corresponds to the pole of free propagators ko = ek — like in the scalar case. The figure [TOl 
(right) represents the | k \ dependence of the amplitudes (pi for a fixed value ko = 0.04. For 
this choice of arguments, the amplitudes are smooth functions of | A; |, though they will be also 
singular for ko > ^ = 0.05. Notice that all the infinities in the BS amplitudes (pi come from 
the free propagator. However the BS amplitude has other non-analytical points. For instance it 
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obtains an imaginary part for it k)"^ > (m + fi)^. 
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6 Electromagnetic form factor 

We demonstrate in this section the advantages of using the BS ampUtude in Minkowski space 
for computing the electromagnetic (em) form factor shown in fig. [TTJ For simplicity, we will 




O(k.P) k o(k,p') 

Figure 11. E.m. vertex in terms of the BS amplitude. 



restrict ourselves to the spinless case for which one has the following expression in terms of the 
Minkowski BS amplitude $m 

(p + pYFM{Q') = -i I ■0^ip + p'-2krirn^-e)^M{l-k,p)<PM(^^-k,p'^ (18) 

As mentioned in the Introduction, the Euclidean BS amplitude cannot provide the right em 
form factor. This fact has been discussed in some detail in our previous work \37\ [38] where two 
kinds of reasons were advocated. 

The first reason was the impossibility for performing the Wick rotation in the integral ()18p . 
Indeed this transformation, represented schematically in figure [T2| consists in replacing the 
integral over the ko variable along the real axis [—R,+R] by the integral along the imaginary 
axis [—iR, +iR]. The result is unchanged provided no any singularity of the integrand is crossed 
when rotating anticlockwise the real into the imaginary integration interval, that is when there 
are no singularities in the first or third quadrants. This is the case, for instance, of those poles 
associated with the free propagators (denoted by + in fig. [T2]) which are located in the second 
and fourth quadrants. It turns out however that the integrand of (|18p has singularities in the 
ko variable (denoted by X in fig. [T2]l which, for some values of fc^, are located in the first or 
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kO=ik4 




kO 



-iR 



Figure 12. Wick rotation in the form factor integral. 



third quadrants. That prevents from performing a naive Wick rotation, i.e. by simply replacing 
in equation (fT8]) ko — )• ik^ and <Pm — ^ "^E, the Euclidean BS amplitudes. Taking properly into 
account the contribution of these additional singularities in the contour integral would require a 
careful analysis on the complex plane of the form factor integrand, a task which is hardly doable 
even by assuming that $m and are known numerically. 

This problem is already manifested at zero momentum transfer in the simplest model case 
of a BS amplitude given by the product of free propagators corresponding to the external lines: 

<Pm = (19) 

[A;2 — + ie] [{p — ky — m? + ie] 

In this case, the integral (jlSp corresponds to the Feynman amplitude represented in fig. Illlwith 
constant vertices [g = 1). The form factor is obtained by multiplying its result (|18p by {p + p')^^ . 
For p = p' it reads: 

AM Fm{0) =zJ (2^)4 lk^_m2 + ie][ip-ky-m^ + ^e]^ ^^^^ 

Since the BS amplitude (fT9|) ~ and consequently Fm - are not normalized, eq. pOj) determines 
the normalization factor. Using the Feynman parametrization 

1 _ /"^ 2xdx 

Jq [(1 - x)a + xhf 

the integral (j'20|) can be calculated analytically and reads: 

1 / Am? arctan , ^ , \ 

which for m = 2 and M = 3 takes the numerical value Fm{0) = 4.99241 x 10~^. 
In the rest frame p = the form factor (j20p has the form: 

^^'^^^ ^ j2n)Hl I ^vr/c^d/c J^^ _ ^2 _ ^2 + [(^^ _°]v//)2 _ ^2 _ ^2 ^ ^,]2 (22) 

In order to check the validity of the " naive Wick rotation" in the form factor integral (|22p , 
we replace there fco = ^^4 and integrate over /c4 in real infinite limits. We obtain in this way: 

-i 2 (M - iki)idki 

- (2^ io i-oo [/fc| + fc2 + ^2]p^ + ^M)2 + fe2 + ^2]2 (23) 
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We have omitted ie in denominator, since the integrand in (j23p is now non-singular at real /C4. 
For particular values of the parameters this integral can be computed numerically. For m = 2, 
M = 3 we get: F/v^(0) = 3.15404 x 10"^ which differs from the value Fm{0) obtained above, 
using Minkowski BS amplitudes. This shows, by a simple example, that performing the Wick 
rotation to the variable /cq - the argument of the BS amplitude ~ in the expression for the em 
form factor is not allowed. 

The reason for this difference is the existence of a pole in the integrand of (j22p created by 
the zeros of the the denominator rightest factor which are given by ko = Mt Vfc^ +m2±ie. If 

< — m^, the pole corresponding to 

ko = M - ^/W+~rn? + ie, 

lies in the first quadrant and it is crossed by the Wick rotation. The residue at this pole is 

-i{M - Vk^ + m2) 

Res = I = 

257r2MVA;2 + m2(2Vfc2 + - M)2 

and its contribution to the form factor is given by: 

FRes{0)= {2iTi Res) Ank'^dk (24) 

Jo 

Numerical calculation for m = 2 and M = 3 gives FRes{0) = 1.83837 X 10"^ and the sum of 
both contributions is in perfect agreement with the value -Fjvf(O). 

Fnw{0) + FResiO) = 3.15404 • 10"^ + 1.83837 • 10"^ = 4.99241 • 10""^ 

This example illustrates, hopefully clearly, how the singularities in the first quadrant indeed 
appear and that they prevent from the Wick rotation. Taking properly into account their con- 
tribution restores the form factor value. However, this cannot be done if one knows only the 
numerical values of the Euclidean BS amplitude. 

The second reason results from the fact that even neglecting the above explained impossibility 
of the Wick rotation in the variable ko, i.e., neglecting the contribution of the singularities 
denoted by X in fig. [121 we still cannot compute the form factor in terms of the Euclidean BS 
amplitude at rest. Indeed, the amplitude ^{k,p) = <P{k'^,p-k) depends on the scalars k'^ = k^ — k^ 
and p-k = poko — pk. After Wick rotation k^ = ik^, the first scalar becomes k'^ = — (A:| -|- k'^). 
Suppose that we are working at the rest frame p = 0. Then the second scalar turns into p-k = 
Poko = iMk4. This purely imaginary value is just the argument of the Euclidean BS amplitude 
and we do not cause any problem. 

However, for non-zero momentum transfer = —{p — p')'^ we have p' 7^ 0. The argument 
p'-k of the amplitude ^{k,p') - right vertex in fig [H] - transforms in the Wick rotation as 
p'-k = p'^ko —p' k = ip'ok/j^ —p' k, i.e. a complex value. The amplitude with complex - not purely 
imaginary - value of its argument p'-k is not the Euclidean one, which enters the left vertex. 
Therefore, even after performing the invalid Wick rotation, the form factor integrand, cannot 
be expressed via BS amplitude obtained by real boost from the Euclidean BS amplitude at rest. 

To avoid the complex boost one can solve the Euclidean BS equation in a moving frame 
with p ^ and obtain in this way the BS amplitude for any value of p. This approach was 
followed in [391 [50] in the framework of the quark model. The form factor was computed in 
terms of the Euclidean amplitudes at p 7^ but without taking into account the contribution of 
the singularities discussed above. 
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Another approximate way to compute the form factors, widely used in the hterature, is the 
so cahed static approximation |36j . It consists in replacing the BS amplitude in a moving frame 
^{k;p'Q,p') by the rest frame amplitude <P{k;pQ = M,p' = 0) boosting only the spatial part of its 
argument k, i.e. letting unchanged the time-component /cq. This means in practice that the left 
vertex of fig. [11] involves the Euclidean BS amplitude withp = 0, ^Eik^, k), while the right vertex 
involves <pE{ki,k + p')- In this way, the form factor is approximately expressed via Euclidean 
BS amplitude. The accuracy of static approximation was estimated in [36] perturbatively, at 
relatively small . 

When using the Minkowski BS amplitude, represented via Nakanishi integral ([2]), we can 
calculate the integral of d^/c analytically and express the form factor in terms of the computed 
Nakanishi weight function (7(7,2). The procedure is now theoretically safe and the precision of 
the result depends only on the accuracy of the numerical solution of the 5(7, z) functions. 

The exact results [371 I38j . together with the static approximation, are shown in fig. [13] 
(left). The difference between the exact calculation and static approximation is small at small 
momentum transfer but it strongly increases with increase of . 

In fig. [13] (right) we compare the exact (Minkowski) form factor with the one found through 
the LF wave function At ah they are almost indistinguishable from each other. 



Euclidian space (sta 


ic approximatLoii) 
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Minkowski space calculation - cross-ladder 
LFD calculation - cross -ladder 



(Q/m)- 



10 20 30 40 50 60 70 80 90 100 
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Figure 13. Left: Form factor via Minkowski BS amplitude (solid curve) and in static approximation (dashed). 
Right; Form factor via Minkowski BS amplitude (solid curve, the same as at left panel) and in LFD (dot-dashed). 



7 Conclusions 

We presented a new method for obtaining the solutions of the Bethe-Salpeter equation in 
Minkowski space, both for two spinless bosons and for two fermions. It is based on a Nakanishi 
integral representation of the BS amplitude and on subsequent Light-Front projection. 

The binding energies for these systems are calculated. They coincide with the ones found via 
Euclidean space solution, thus providing a validity test of our method. For two fermions, the 
solutions for the scalar exchange and positronium states without vertex form factor [A — )■ oo) 
are stable below some critical value gc of the coupling constant, respectively gc = 27r (scalar 
exchange) and gc = (positronium). 

The BS amplitudes are obtained in terms of the computed Nakanishi weight functions. They 
exhibit a singular behaviour due, on one hand, to the poles of the free propagators and, on 
the other hand, to non analytic branch cuts which are responsible for the appearance of an 
imaginary part. 
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Some applications to em form factors are also presented. We demonstrate that the Wick 
rotation in the variable ko - argument of the BS amplitude - is not allowed. It is prevented 
by the singularities which appear in the first quadrant of the complex /cQ-plane and whose 
contributions are hardly evaluable in practice. Neglecting these contributions, one obtains an 
approximate expression which supposes the knowledge of Euclidean BS amplitude in a non- 
zero momentum frame. The latter requires to solve numerically a much more involved - though 
solvable - equation than at rest. On the contrary, in terms of the BS amplitude in Minkowski 
space, represented in the Nakanishi form ([2]), the form factors can be calculated without any 
problem. 



Appendix A: Numerical methods 



The numerical methods wil be illustrated by considering the general case of Ua coupled two-dimensional integral 
equations for the Nakanishi weight functions ga in the form: 

/ d7' / dz' V^{'y,z,'y',z') ga{'y',z') = ^ H f dz' Vaa'{l,z,^' ,z') Qa'i^i' ,z') (A.l) 

Jo J -1 ^ ^1^-^ Jo J -1 

For the scalar case the number of amplitudes is = 1 while for fermions it depends on the quantum numbers of 
the state (na=4 for J = 0, Ua = 8 for J — 1). 

Equation (jA.l[) has been solved by expanding the unknown functions Qa on a bicubic spline basis [21] over a 
compact integration domain O — x I~ = [0,^max] ^ [^^i+l]- 



2N^ + l 2N^+1 
i=0 j=0 



(A.2) 



The interval corresponding to variable x = 'y,z is divided into subintervals - [a;i, a;i+i] with i = Q, . . . - 
satisfying xo = Xmin and xn^ = Xmax- The distribution of grid points Gx ~ {xi} are adjusted to the structure of 
the solution. By means of expansion (|A.2|) . equation (jA.l|l can be written in the form 

(A.3) 



with 



Ba'i'j'{'y,z) = 



d-y' / dz' V^{'y,z,'y',z') Si'{j')Sf{z') 5aa' 



d^' I dz' Vaa'{'y,Z,^',z') Si^{'j')Sj^{z') 



(A.4) 
(A.5) 



The left hand side of (|X3| is in fact diagonal in the number of amplitudes. 

The spline functions for the expansion on a variable x, Sj{x), are defined once the corresponding grid points 
Gx = {xo,xi, . . . , xjVj: } are fixed. These are 2Nx + 1 functions, usually denoted hy So, Si, . . . S'2]Vj;+i. 

The functions S2i and 52^+1 are associated to the grid point Xi and have a support limited to the two 
consecutive intervals surrounding it, i.e. [xi-i, Xi+i], a property that makes easier the computation of integrals 
(|A.4[) and (|A.5[l . Their analytic expressions are given by 



S2i{x) 



S2t+l{x) 



+ 



+ 



X — Xi 




Xl — X-i 












x — x^ 




Xl—X 






=&) 











if a; G [x^, Xi+i] 
if a; ^ [xt-i,x,+ i] 

{xi — Xi-i) X e [xi-i,Xi] 

(xi+i - Xi) X G [xi,Xi+i] 

if X ^ [x,-i,x,+ i] 



and they are such that the values of the solutions at the grid points are simply given by the even coefficients of 
the expansion (|A.2p . i.e. : 

gai'yi,Zj) = ga,2i,2j (A. 6) 
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PIH cubiques 



Figure 14. Cubic spline functions associated to the grid point Xi. 



They are represented in figure 14. 

Equation (|A.3[) is validated on a ensemble {7;, Zj} C O of {2N~/ + 2){2Nz + 2) "collocation" points which are 
taken equal to the N=2 Gauss quadrature abcises inside each subinterval of Ix plus the 2 borders. This leads to 
a generalized eigenvalue problem 

A B(M)g = A{M)g (A.7) 

in which matrices B and A represent respectively the integral operators of the left- and right-hand sides of (|A.1|) 
and have the matrix elements 

^aij,a'i' j' -^a' i' j' 5 -^J ) 

The unknown coefficients Qaij are determined by the solutions of (|A.7|) corresponding to A(A/) = 1. 
Notice that the value of the solutions at the collocation points are given by the matrix product 

<7a(7i,%) ~ ^ ^ 9a'i'j' 
i'j' 

with the spline matrix U given by 

Uij.i'j' = Sii ('yi,)Sj' (zj) 

It turns out that the discretized integral operator B has very small eigenvalues. They are unphysical but make 
unstable the solution of the system (|A.7|) . To regularize B, we have added a small constant er to its diagonal 
part on the form: 

B^B + erU (A.8) 

This procedure allows us to obtain stable eigenvalues with an accuracy of the same order than er until values of 
Er as small as 10^^^. 

The property HA.6I) is very useful to implement boundary conditions to the solutions we are interested in. 
Although not necessary when working in momentum space integral equations, it uses to help the convergence of 
the numerical algorithms and allows to avoid unwanted solutions, which can be either due to the discretisation 
on a finite box or having other kind of symmetries. If we know in advance, for instance, that the solution vanishes 
at 2; = ±1, this implies according to (|A.6|) to impose to the coefficients ga,i,o = ga,i,2N = 0, that is in practice to 
remove them from the expansion (|A.2|) with the consequent reduction in the dimension of the matrix equation 
(Ell- 
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